Relativistic BCS-BEC Crossover at Finite Temperature 
and Its Application to Color Superconductivity 
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The non-relativistic GoG formalism of BCS-BEC crossover at finite temperature is extended 
to relativistic fermion systems. The uncondensed pairs contribute a pseudogap to the fermion 
excitations. The theory recovers the BCS mean field approximation at zero temperature and the 
non-relativistic results in a proper limit. For massive fermions, when the coupling strength increases, 
there exist two crossovers from the weak coupling BCS superfluid to the non-relativistic BEC state 
and then to the relativistic BEC state. For color superconductivity at moderate baryon density, the 
matter is in the BCS-BEC crossover region, and the behavior of the pseudogap is quite similar to 
that found in high temperature superconductors. 
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I. INTRODUCTION 

It is well-known that, by adjusting the attractive cou- 
pling strength among the constituents, a fermion sys- 
tem may undergo a smooth crossover from the Bardeen- 
Cooper-Shriffcr (BCS) superfluidity/superconductivity 
in degenerate fermion gas to the Bose-Einstein conden- 
sation (BEC) of composite molecules. Such a BCS-BEC 
crossover is theoretically due to the fact that the wave 
functions of BCS and BEC ground states are essentially 
the same[l|,0]- The BCS-BEC crossover is expected to be 
realized in high temperature superconductor and atomic 
fermion gas [JH, H, [g, 0, @| via using an external magnetic 
field to change the s-wave scattering length Q. 

The superconductivity in Quantum Chromodynamics 
(QCD), i.e., the color superconductivity fiol] . is naturally 
considered as a system to study the relativistic BCS- 
BEC crossover. Due to the asymptotic property of QCD, 
there may exist a crossover from the BCS superconduc- 
tivity with weakly bound quark pairs at high baryon 
density to the BEC state of composite hadrons at low 
baryon density fill]. Such a BCS-BEC crossover in QCD 
may also be realized in chiral condensed matter [l^, EH 
and in pion superfluid [TTI] . At moderate baryon density, 
while a diquark BEC state may not be realized due to 
the chiral symmetry restoration, the attractive coupling 
strength is obviously not located in the weak coupling re- 
gion. It is shown in many effective QCD models that, the 
quark en ergy gap at moderate baryon density is about 
100 MeVpjl] which is already of the order of the Fermi 
energy. The strong coupling in this case may induce a 
so-called pseudogap effect, which has been investigated 
in two flavor color superconductivity above the critical 
temperature pl| . A natural question is how the pseudo- 
gap modifies the critical temperature and thermodynam- 
ics of the color superconductor. To answer this question, 
one needs to construct a relativistic theory at finite tem- 
perature which can describe the pseudogap and possible 
BCS-BEC crossover. 

The BCS-BEC crossover in relativistic fermion systems 
are recently investigated in the Nozieres-Schmitt-Rink 



(NSR) theory above the critical temperature [r| [13], 
the boson-fermion model[H| and the BCS-Leggett mean 
field theory at zero temperature[19]. It is shown that, 
not only the BCS superfluidity and the non-relativistic 
BEC (NBEC) of heavy molecules but also the NBEC 
and the relativistic BEC (RBEC) of nearly massless 
molecules can be smoothly connected. In the RBEC 
state, anti-fermion pairs (anti-bosons) are excited and 
become nearly degenerate with fermion pairs (bosons). 
From the NSR theory at T > T c , where T c is the critical 
temperature, the difference between the NBEC[2lj and 
RBEC [H H| states is significant [U, E3 • 

It is widely known that, at zero temperature the mean 
field theory is a good approximation to describe the BCS- 
BEC crossover |2fjj. and the pair fluctuations can be safely 
neglected even at strong coupling. Only around the uni- 
tary limit, i.e., the infinite scattering length limit, the 
pair fluctuations are somewhat important to obtain a 
proper value of the universal constant 0- In our pre- 
vious paper we investigated the generalization from 
non-relativistic to relativistic BCS-BEC crossover at zero 
temperature in the BCS-Leggett mean field theory. At fi- 
nite temperature, however, the condensed pairs with zero 
momentum can be thermally excited, and one should go 
beyond the mean field approximation to treat properly 
the uncondensed pairs [6]. 

There exist many methods to treat pair fluctuations 
at finite temperature. In the NSR theory, which is also 
called GqGo theory, the pair fluctuations enter only the 
number equation, and the fermion loops which appear in 
the pair propagator are constructed by bare Green func- 
tion Go- As a consequence, such a theory is in principle 
not self-consistent and is valid only at T > T c . For the 
study of BCS-BEC crossover, one needs a theory which is 
valid not only above the critical temperature but also in 
the symmetry breaking phase. While such a strict theory 
is not yet reached so far, some T-matrix approaches are 
recently developed, see for instance [f| [2(| ■ Among them, 
the asymmetric pair approximation or the so-called GoG 
scheme [(I l24j is a competitive one. The effect of the pair 
fluctuations in the GqG method is treated as a fermion 
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pesudogap which has been widely discussed in high tem- 
perature superconductivity. In contrast to the NSR the- 
ory (Go Go scheme), the GoG scheme is self-consistent 
and keeps the Ward identity Q. 

In the study of color superconductivity at moderate 
density, the color condensed phase is of great interest. 
The NSR theory [l^, [It}, which seems valid in the nor- 
mal phase, can only predict the transition temperature of 
color superconductivity. A necessary task in this field of 
research is to develop a relativistic BCS-BEC crossover 
theory in the symmetry breaking phase. In this paper, 
we will generalize the GoG scheme to relativistic fermion 
systems. A necessary requirement for such a generaliza- 
tion is to recover the non-relativistic limit 0] and mean 
field limit Q properly. With this theory, we can calculate 
the critical temperature T c for arbitrary coupling and 
describe the BCS-NBEC-RBEC crossover at finite tem- 
perature. It, as an application, can be used to study the 
pseudogap effect on color superconductivity. 

The paper is organized as follows. In section |TT] we 
review the BCS mean field theory for relativistic super- 
fluidity/superconductivity. In the framework of the GoG 
scheme, we include in Section IIIII the contribution from 
the uncondensed pairs and construct coupled equations 
for the superfluid order parameter and pseudogap. In 
section IIV1 we apply the theory to massive fermions and 
study the BCS-NBEC-RBEC crossover at finite temper- 
ature. In section [Vj we apply the theory to color su- 
perconducting quark matter. We will calculate the tran- 
sition temperature and the quark pseudogap and show 
the significance of the fluctuations at moderate baryon 
density. We summarize in section IVII 



II. BCS MEAN FIELD THEORY 

We consider a model with only fermions as elementary 
blocks. The Lagrangian density can be written as 



C = tp (i^d^ — m)ip + Ci 



(1) 



where ip,ip denote the Dirac fermion fields with mass 
77i, and Cj indicates the attractive interaction among 
fermions. Since the dominant interaction is the J p = + 
scalar channel, the interaction for the pairing between 
different spins can take the form[l(| 



(2) 



where g is the attractive coupling constant, and G = 
*7o72 is the charge conjugation matrix. Generally, by ad- 
justing the coupling strength, the crossover from conden- 
sation of spin-zero Cooper pairs with large size at weak 
coupling to the Bose-Einstein condensation of deeply 
bound bosons at strong coupling can be realized. In 
our model, only fermions are elementary particles. An- 
other type of model which is used to discuss the BCS- 
BEC crossover in high temperature superconductors and 



atomic Fermi gases is the so-called boson-fermion model 
where both fermions and bosons are considered as ele- 
mentary blocks. Such a model is recently generalized to 
study the relativistic BCS-BEC crossover |18j. 

In order to develop a finite temperature theory includ- 
ing pair fluctuations in the symmetry breaking phase, we 
first review in this section the BCS mean field theory in 
the functional integral approach and GqG formalism. 



A. Functional Integral Approach 

In the functional integral approach, we start the cal- 
culation from the partition function in imaginary time 
formalism, 



d 3 x(£+Ki/> t 'iA) 



(3) 



where j3 is the inverse temperature, j3 = 1/T, and 
(i is the chemical potential corresponding to the net 
charge density ip^ip and is determined by the charge 
conservation. Performing a Hubbard-Stratonovich trans- 
formation which introduces an auxiliary pair field 
A(x) = g4> T ' (x)Ci r y5t(j(x) /2, and then integrating out the 
fermions, we derive the partition function 



Z = / DADA* e - s <= f f[A,A*] 



with the effective boson action 



S e fi = 



P 



dT 



d 3 , 



9 2/3 



(4) 



(5) 



in terms of the inverse Nambu-Gorkov propagator 

G" 1 = ij^dfj, - ?77 + fjrf a 3 + 775A0 + + 775 AV_, (6) 

where a± = [a\ ± ia<x)j1 are defined in the Nambu- 
Gorkov space with <7j(i — 1,2,3) being the Pauli ma- 
trices. 

The mean field theory is a good approximation to de- 
scribe the BCS-BEC crossover at low enough tempera- 
ture, namely T <C T c , since the dominant contribution of 
fluctuations to the effective potential is from the Gold- 
stone mode and is proportional to T 4 Q. In the mean 
field approximation, we consider a uniform static saddle 
point A (a;) = A sc which satisfies the stationary condi- 
tion SS c s[A sc ]/5A sc = 0. The thermodynamic potential 
fimf = S c g[A sc )/(pV) at the saddle point can be evalu- 
ated as 



fin 



d 3 k r 



(2tt) 3 



ln(l 



-0K 



ln(l 



-0E. 



>) 



(7) 



where we have defined the quasi-particle energies E^ = 



±)2 + A2 c with e 



£k ± M and ek = Vk 2 + m 2 . 
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Minimizing S7 m f , we get the gap equation to determine 
the order parameter A sc in the symmetry breaking phase, 



1 
9 



d 3 k 

(2^)3 



l-2/(£ k ) , 1-2/0E+) 



2E,: 



2Ei 



(8) 



where f(x) = \j{eP x + 1) is the Fermi-Dirac distribution 
function. In the study of BCS-BEC crossover, people 
often consider the thermodynamics in canonical ensem- 
ble with fixed fermion density n by fixing the Fermi mo- 
mentum kf through the relation n = kf/(3ir 2 ) at zero 
temperature. At finite temperature, the density can be 
obtained from the first order derivative of the thermody- 
namic potential with respect to the chemical potential, 



d 3 k 

(2tt) 3 



l-i^(l-2/(££)) 



1-^(1-2/(0) 



(9) 



The first and second terms in the square bracket on the 
right hand side of equations (HJ) and ([9]) correspond re- 
spectively to fermion and anti-fcrmion degrees of free- 
dom. 



B. GqG Formalism 

Now we reexpress the BCS mean field theory in the 
GqG formalism 0, 0, HE]- Such a formalism is convenient 
for us to go beyond the BCS and include uncondensed 
pairs at finite temperature. Let us start from the fermion 
propagator S in the symmetry breaking phase. The in- 
verse propagator reads 



with the inverse free propagator 

QQ X {k, = (iuj n + ^)7o - 7 • k - m, 



(10) 



(11) 



where k = (zc<j n ,k) is the fermion four momentum at 
finite temperature with uj n being the fermion frequency 
Lu n = (2n + 1)ttT (n = 0, ±1, ±2, ■ • • )■ The propagator 
can be formally expressed as 



S(k) = 



with the diagonal and off-diagonal elements 

g(k,n) = [g - 1 (fc,^)-s sc (A ; )]" 1 , 

T{k,^i) = -g(k,fi)i'y 5 A sc go(k,-^), 

where the fermion self-energy S sc is defined as 

£ sc (fc) = ij 5 A sc g (k,-^)ij 5 A sc 
= -A 2 sc g (-k,»). 



(12) 



(13) 



(14) 



With the help of the energy projectors 



A±(k) 



1 



1± 



70 (7 • k + to) 



(15) 



the propagator elements can be explicitly evaluated as 

(*^n + C k ) A +70 (^n-Ck)A_7o 



g(k, n) = 



F{k,n) 



iA sc A+7 5 



«A SC A_7 5 



(ico n r-(E k r ^„) 2 -«)= 



(16) 



The gap equation for the order parameter A sc is related 
to the off-diagonal element, 



■ fJ 



|A sc ^Tr[S(fc, m)], (17) 

k 

and the fermion number is controlled by the diagonal 
element, 



= -i^2^[-y g(k,fi)} 



(18) 



with the four momentum integration J^k = 
iTJ2 n I ^ 3 k/ (27r) 3 at finite temperature. Com- 
pleting the Matsubara frequency summation, we can 
reobtain the gap equation (|8|) and number equation (|9j). 

In the BCS mean field theory, fermion-fermion pairs 
and anti-fermion-anti-fermion pairs explicitly enter the 
system below T c only through the condensate A sc . In the 
GqG formalism, the fermion self-energy can equivalently 
be expressed as 



s sc(fc) = y^4c(g)0o(g - k,fi) 



(19) 



associated with a condensed-pair propagator given by 

(20) 



where q — (iv n ,q) is the boson four momentum with 
boson frequency v n — 2nnT. 

The BCS theory can be related to a specific pair sus- 
ceptibility x defined by 



Xbgs(q) = -^^2Tr{g(k,fi)g (q- k,fi)} 



(21) 



with which, the gap equation for the condensate A sc can 
be written as 



I-SXbcs(O) = 0. 



(22) 



This implies that the uncondensed pair propagator 
should be of the form 



t{q) 



i<J 



1 -5XBCs(g)' 



(23) 
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and t~ 1 (q — 0) is proportional to the pair chemical po- 
tential /ipair. Therefore, the fact that in the symmetry 
breaking phase the pair chemical potential is zero leads 
to the BEC-like condition 

t~ 1 (q = 0)=0. (24) 

While the uncondensed pairs do not play any real role 
in the BCS mean field theory, such a specific choice of 
the pair susceptibility and the BEC-like condition tell us 
a way how to go beyond the BCS mean field theory and 
include the effect of uncondensed pairs. 

III. BEYOND MEAN FIELD THEORY 




FIG. 1: Diagramatic representation of the propagator t pg for 
the uncondensed pairs and the fermion self-energy. The to- 
tal fermion self-energy contains contributions from condensed 
(Esc) and uncondensed (S pg ) pairs. The dashed, thin solid 
and thick solid lines in £ pg represent, respectively, the cou- 
pling constant g/2, bare propagator Go and full propagator 
Q. This diagram is taken from [291 ]. 



While the uncondensed pairs can be safely neglected 
at weak coupling, they should be included for a self- 
consistent theory at arbitrary coupling and at finite tem- 
perature. We now go beyond the BCS mean field ap- 
proximation and include the uncondensed pairs in the 
GqG formalism. It is clear that, in the BCS mean field 
approximation the fermion self-energy S sc includes con- 
tribution only from the condensed pairs. At finite tem- 
perature, the condensed pairs with zero total momen- 
tum can be thermally excited, and the total propagator 
should contain both the condensed (sc) and uncondensed 
or "pseudogap" -associated (pg) contributions, 



t(q) = t ps (q) +t sc (q), 



W<?) 



If) 



1 - 9x(q) ' 



t sc (q) = i^S(q). 



(25) 



Now the total fermion self-energy becomes 



E(fc) 



t(q)Go(q ~ k, m) = S sc (fc) + E pg (fc), (26) 



with the mean field part 

£ BC (fe) =J2t sc (q)g (q-k,fj,) (27) 
i 

and the pseudogap related part 

Epg(*)=X>g(9)0o(g-*,/i). (28) 
i 

With the full propagator 

g(k, f i) = [g 1 (k,(i)-x(k)]- 1 (29) 

in terms of the total self-energy, the pair susceptibility is 
still given by 



x(q) 



■|^Tr^(fc )M )eo(9-fc,M)- (30) 

k 



The GqG formalism used here is diagrammatically illus- 
trated in FigfT] 

Note that, the feedback of the pair fluctuations on 
the order parameter A sc is included, and it and the 
chemical potential /j, are in principle determined by the 
BEC condition t p ~ g 1 (0) = and the number equation 

n = -iEfc Tr [7o£(&,Ai)]- 

The above equations are hard to handle analytically. 
In the symmetry breaking phase with T < T c , the BEC 
condition <pj(0) = implies that t pg (q) is peaked at 
q = 0. This allows us to approximate 



E(A;)~-A 2 So(-/c,M), 



(31) 



where A 2 contains contributions from the condensed and 
uncondensed pairs, 



A = Af 



pg 



with the pseudogap A pg defined as 



E 



*pg(?)- 



(32) 



(33) 



It is necessary to point out that, above the critical tem- 
perature T c such an approximation is no longer good, 
since the BEC condition is not valid in normal phase. 

Since the pseudogap A pg looks similar to the conden- 
sate A sc , a natural question is whether a finite A pg breaks 
the symmetry of the system. If yes, A sc will no longer 
be considered as the order parameter of the phase tran- 
sition. By omitting a term of the order of O (A 2 C /A 2 ), 
where A is a momentum cutoff, the inverse fermion prop- 
agator including the feedback of the pair fluctuations can 
be written as 



-i (k \ _ (Go \k,fj,) - S pg (fc) 

1 ] \ *7 5 a sc g -\k, 



where £' 



«75A SC 

-M) - K g (k) 

(34) 

Jp ^pg(^ - * — m)- It ls now clear that, 
the pseudogap appears in the diagonal elements of the 
Nambu-Gorkov propagator and does not break the sym- 
metry of the system. On the other hand, parallel to the 
discussion in non-relativistic theory [6ll7l.l24l]. we can show 
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that Ap g is just the fluctuation of the order parameter 
field A(x), 

|A|) 2 , (35) 



A 2 g = (|A| 2 )- 



and hence it does not break the symmetry. 

Under the approximation pip , all the equations in 
the mean field theory are still valid, the only change 



is the replacement of E^ = y + A 2 C by E± = 

J (£k) 2 + A 2 . For instance, the diagonal element Q of 
the full propagator, the fermion number n and the gap 
equation for A take exactly their mean field forms ([To]) . 
© and ||SJ). The equations ©, © and (33]) determine 
self-consistently the order parameter A sc , the pseudogap 
Apg and the chemical potential fj, as functions of tem- 
perature T. Note that the pair fluctuation effect is self- 
consistently included in the coupled equations through 
the pseudogap A pg . It is necessary to point out that, the 
GoG approach we used is quite different from the NSR 
theory. In the NSR theory, the pair fluctuations enter 
only the number equation via adding a molecule number 
termpE[l7t. 

However, solving such a coupled set of equations is 
still rather complicated. Fortunately, the BEC condition 
allows us to do further approximations for the pair prop- 
agator t pg (q). Using the BEC condition 1 — gx(0) = 0, 
the T-matrix can be written as 



*pg(«) 



x(q) - x(o) 



(36) 



Since the pseudogap is dominated by the gapless pair 
dispersion in long wavelength limit, we can expand the 
susceptibility around q — in this limit, 



W<?) 



Z iqa + Z 2 ql-eci 2 +iT{qy 
where the coefficients Z\ , Z 2 and £ 2 are defined as 



(37) 



Z = — 
1 dq 



Zi = - 



ld 2 X 



q=0 
2 



t l _ _ " A. 



2dq 2 



q=0 



(38) 



and we have considered the fact that the susceptibility 
depends only on q 2 . The explicit expressions for Z\, Z 2 
and £ 2 are listed in Appendix |A"1 

In the symmetry breaking phase where the tempera- 
ture is low, it is believed that the pairs are long-lived 
and we can neglect their width T. With the expansion 
for the pair propagator, the equation (|33[) now takes a 
simple form 



1 

Z~2 



d 3 q 1 + b(ujq -v) + i)(w q 
(2tt) 3 2uj^ 



(39) 



where b(x) — \ / (el 3x ~X) is the Bose-Einstein distribution 



function and cj q and v arc defined as 



c 2 q 2 , 



2Z 2 



Z 2 



(40) 



The first term on the righthand side of (|39f suffers 
ultraviolet divergence, but it can be dropped out via 
renormalizationjy, 0] ■ 

Let us first discuss some conclusions from the above 
equations without detailed numerical calculations. 

1) At zero temperature, the pseudogap A pg vanishes au- 
tomatically and the theory is reduced to the BCS mean 
field approximation 19] . 

2) If the coupling is not so strong that the molecule bind- 
ing energy E b satisfies E b <C 2m, the theory is reduced 
to its non-relativistic version Q for systems with kf <C m 
or n <C m 3 . 

3) If Z\ dominates the propagator i pg , the pair dispersion 
is quadratic in |q|, and the pseudogap A pg can be an- 
alytically integrated out and is proportional to T 3 / 4 at 
low temperature. On the other hand, if Z 2 is the dom- 
inant term, the pair dispersion is linear in |q| and A pg 
becomes proportional to T at low temperature. In the 
next section, we will show that the first case happens 
in the NBEC region and the second case occurs in the 
RBEC region. 

4) From the explicit expression of Z\ shown in Appendix 

m 



n f d 3 k 
2 ~ J (27p 

brackets is just the total num- 

nrl tiqii-g fKncrtnel 



(/(C) (41) 



the quantity in the square brackets is just tJ 
ber density of the bound pairs (bosons) 



"B 



= ^A 2 



(42) 



From the relation A 2 = A 2 C + A pg , ub can be decom- 
posed into the condensed pair number n sc and the un- 
condensed pair number n pg , 

?i sc = ZiA 2 c , n pg = Z!A 2 pg . (43) 

The fraction of the condensed pairs can be defined by 

A 2 



Pc = 



n/2 



(44) 



5) In the weak coupling BCS region, we expect the 
fermion number density 



d 3 k 

(2T) 3 



(/«k)-/(^)) 



(45) 



which leads to = in this region. In the deep BEC 
region, however, almost all the fermions form two body 
bound states which results in 7Jb — n/2. At zero temper- 
ature, we have A pg = 0, tib = n sc and P c = 1, while at 
the critical temperature T c , the order parameter A sc dis- 
appears, and the uncondensed pair number n pg becomes 
dominant and is approximately equal to n/2. 

Numerically, the transition temperature T c can be cal- 
culated from (|3"9"|) and the generalized equations ([5]) and 
© by setting A sc = 0. Usually, at and above T c where 
the order parameter A sc disappears, the pseudogap A pg 
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does not vanish. We can define a limit temperature T* 
where the pseudogap starts to disappear. Between the 
two temperatures T c and T* is the so-called pseudogap 
phase. While the present G$G formalism is likely valid 
only in the symmetry breaking phase with T < T c , it can 
be generalized to the region above T c by introducing a 
non- vanishing pair chemical potential /x pa i r [3|- We will 
do such a generalization, but the numerical results in the 
following will be presented mainly at T < T c . 

Above the critical temperature T Cl the order parameter 
A sc vanishes, and the BEC condition is no longer valid, 
1 — 5x(0) 7^ 0- As a consequence, the propagator of the 
pair takes the form 



X (q) - X(0) - Z Q 
x(0). As an estimation of A t 



(46) 



and 



With Zq = l/g - J^a c.i oniiuouiuu ui ^-> p g 

T* , we still perform the expansion for the susceptibility 
around q = 0, 



Z iqo + Z 2 q%-e\q\ 2 -Zo + ir(q)- 
Now the pseudogap equation becomes 



A z = — 

pg 



1 

Z~2 



(2tt) s 



2^ q 



with the definition 



(47) 



(48) 



(49) 



The equation (|48|) together with the number equation ([9]) 
determines the pseudogap A pg and chemical potential /i 
above T c . Since the pair dispersion is now no longer gap- 
less in the long- wavelength limit, and Zq will generally 
increase with temperature, we expect that A pg will drop 
down and approach zero at the dissociation temperature 

In the end of this section, we discuss the thermody- 
namics of the system. The naive BCS mean field theory 
does not include the contribution from the uncondensed 
bosons which, however, dominate the thermodynamics at 
strong coupling. In the present theory, considering the 
uncondensed pairs, the total thermodynamic potential 17 
contains both the fermion and boson contributions, 



^ ^cond ~t~ ^fermion ^bosoni 

where S7 C ond is from the condensed pairs, 

A 2 

o — sc 

^ ^cond i 

9 

^fermion from the fermion excitations, 
d 3 k 



n 



fermion 



(2tt) 3 



- /3E k) + l n (l + e 



(50) 



(51) 



(52) 



-0E- 



and Oboson from the uncondensed pairs, 

^boson = ^2 M 1 - 9X(q)}- (53) 

9 

Under the approximation (|37p for the pair propagator, 
the boson part in the symmetry breaking phase can be 
evaluated as 



d 3 q 

(2^)3 



ln(l 



In (1 



with OJq 



(54) 



There exist two limiting cases for the boson contri- 
bution. If Z\ dominates the pair propagator, the pair 
dispersion is quadratic in q, and Jlboson recovers the ther- 
modynamic potential of a non-relativistic boson gas, 



boson 



i 



/3 J (2tt) 3 



In 1 



-/3q 2 /(2m B ) 



(55) 



On the other hand, if Z 2 dominates the pair properties, 
the pair dispersion is linear in |q| and we obtain the ther- 
modynamic potential for an ultra-relativistic boson gas 



9. 



UR 
boson 



rf 3 q 

13 j (2^y- 



ln 



(l-e-W^l) 



(56) 



As we will show below, the former and the latter corre- 
spond to the NBEC and RBEC region, respectively. The 
bosons and fermions behave very differently in thermo- 
dynamics. As is well known, the specific heat C of an 
ideal boson gas is proportional to T a with a — 3/2 and 
3 corresponding to non-relativistic and ultra-relativistic 
systems, but the naive BCS mean field theory predicts 
an exponential law C cx e~ A °/ T , where Ao is the gap at 
zero temperature, A = A(T = 0). 



IV. BCS-NBEC-RBEC CROSSOVER WITH 
MASSIVE FERMIONS 

In this section, we study the BCS-BEC crossover when 
the coupling constant g increases. Since our model is non- 
renormalizable, a proper regularization is needed. In the 
case with massive fermions we employ the often used non- 
relativistic regularization, namely, to replace the bare 
coupling g by a renormalized coupling £/(l7l. [l9|. 



1 

U 



|k|<A 



d 3 k 

(2V) 3 



(57) 



The effective s-wave scattering length a s can be related 
to U by U = A-Ka s /m. While this is a natural extension 
of the non-relativistic regularization to relativistic sys- 
tems, the ultraviolet divergence can not be completely 
removed, and a cutoff A still exists in the theory. In this 
regularization, the solution of the coupled equations de- 
pends on three dimensionless parameters: the effective 
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coupling constant r\ — l/(kfa s ), the quantity £ = kf/m 
which reflects the fermion number density, and the cutoff 
A/to. 

We assume in this section that the fermion density n is 
not very high and satisfies the relation n < to 3 or £ < 1. 
In this case the system is not ultra-relativistic and can 
even be treated non-relativistically in some parameter re- 
gion. From the study in NSR theory above T c and in the 
BCS-Leggett theory at T = 0, if the dimensionless cou- 
pling rj varies from — oo to +00, the system will undergo 
two crossovers fl(I IT7I \l§L . the crossover from the BCS 
state to the NBEC state around rj ~ and the crossover 
from the NBEC state to the RBEC state around 77 ~ £ _1 . 
The NBEC state and the RBEC state are characterized 
by the molecule binding energy Eb. We have Eb <C 2to 
in the NBEC state and E b ~ 2m in the RBEC state. 
1)BCS region. In weak coupling BCS region, there exist 
no bound pairs in the system. In this case, Z\ is suffi- 
ciently small and Z2 dominates the pair dispersion and 
we have A 2 g ex l/(Z 2 c 3 ) after a simple algebra. Since 
A should be small in the weak coupling region, and c 
can be proven to be approximately equal to the Fermi 
velocity 7], the pseudogap A pg is very small and can be 
safely neglected in this region, as we expected. There- 
fore, the BCS mean field approximation is good enough 
at any temperature, and the critical temperature satisfies 
the well known relation T c ~ 0.57Ao. For example, in the 
non-relativistic limit with kf -C to, the anti-fermion de- 
grees of freedom can be ignored and the pair susceptibil- 
ity recovers its non-relativistic version^, see the result in 
Appendix [XJ The critical temperature can be expressed 
asy 



-<7 



,27,/TT 



(58) 



where 7 is the Euler constant and e/ = fc?/(2m) is the 
Fermi kinetic energy. In this region, even though Z2 dom- 
inates the pair dispersion, we can show that c cx A is 
vanishingly small due to the weak coupling. Since the 
boson contribution to thermodynamics can be neglected, 
the specific heat at low temperature takes the well known 
form C cx e - A °/ T . 

2)NBEC region. In the non-relativistic BEC region 
with 77 > 1 but ?] <C C -1 [Ill> the molecule binding en- 
ergy Eb is much less than 2to, namely — m\ <C to, the 
boson mass is approximately 2m, and the system can 
be regarded as a non-relativistical boson gas, if kf/m is 
small enough. Assuming kf -C to, the anti- fermion de- 
grees of freedom can be neglected, and we can recover 
the non-relativistic result [6j. In this region, the gap A 
becomes as large as the Fermi kenetic energy e/. From 
Z\ oc 1/A 2 and Z2 oc 1/A 4 , Z\ is the dominant one and 
the pair dispersion becomes quadratic in |q|. In this case, 
the propagator of the uncondensed pairs can be approx- 
imated by 



iz- 1 



90 



|q|7(2m B )' 



(59) 



where the pair mass TOb is defined by tob = Zi/2£ 2 , and 
we have the simple relation 



m B T\ 



3/2 



2tt 



c 



(60) 

Since ZiA 2 g is equal to the total boson density n B at 
T = T c , we arrive at the standard critical temperature 
for Bose-Einstein condensation in non-relativistic boson 
gas [2l| , 



2tt 
m B 



C(|) 



2/3 



(61) 



The boson mass tob is generally expected to be equal 
to the boson chemical potential u B = 2/i. In the non- 
relativistic limit kf <C to, we can show tob — 2to and 
ZiAp g ~ n/2 at T = T c , the critical temperature be- 
comes T c — 0.218e/. Since Z\ dominates the pair dis- 
persion, the pseudogap is proportional to T 3 / 4 and the 
specific heat is proportional to T 3 / 2 at low temperature. 
3)RBEC region. In this region we have the molecule 
binding energy Eb — * 2to and chemical potential fj, — * 0. 
In this case, non-relativistic limit can not be reached even 
for kf <C to(19|. Since the bosons with mass tob = 2/i 
become nearly massless in this region, the anti-bosons 
and anti-fermions can be easily excited, and the system 
contains both bosons and anti-bosons. From the relation 



n>b ~ n h = ZiA 



(62) 



at T — T c , where rib and are the boson and anti- 
boson numbers, while rib and are both very large, their 
difference produces a small pure boson density tib — n/2. 
On the other hand, for fj, — * we can expand Z\ in powers 
of chemical potential fi, 

Z 1 ~Ru + 0(u 3 ) = ^m B + 0(u 3 ) (63) 

and hence Z2 dominates the pair dispersion, which means 
that the pseudogap is proportional to T at low temper- 
ature. In this case, the propagator of the uncondensed 
pairs can be approximated by 



—iZ 



1 2 o 



c 2 |q| 2: 



which leads to the relation 

d 3 q 6(c|q|) 



Z 2 A 



T 2 



(2tt) 3 c|q| 12c 3 ' 
Combining the above equations, we find 

/24c 3 Z 2 n B N ' 1 



V R 



TOB 



(64) 



(65) 



(66) 



In the RBEC limit /1 — > 0, we can approach to the 
standard critical temperature for ultra-relativistic Bose- 
Einstein condensation 1221 l23i. 



f 3n B 



1/2 



(67) 
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Since hb is almost fixed and me — > 0, T c would approach 
to infinity in the RBEC limit. In the ultra-relativistic 
boson gas, the specific heat at low temperature is pro- 
portional to T 3 . 

We now turn to numerical calculations. From the cou- 
pled equations (JSj> , (J9j) and (|39|) . we can solve the criti- 
cal temperature T c , chemical potential /i(T c ) and pseu- 
dogap A pg (T c ) at T c as functions of the coupling rj at 
fixed kf/m. In Figj2]we show the numerical results with 
the parameters A/to = 10 and kf/m = 0.5. The BCS- 
NBEC-RBEC crossover can be seen directly from the be- 
havior of the chemical potential \i. In the BCS region 
with — oo<f7<0.5,yuis larger than the fermion mass 
and approaches to the Fermi energy Ef in the weak cou- 
pling limit r\ — ► — oo. The NBEC region is located around 
—0.5 < -q < 4 and the RBEC region is at about r] > 4. 
The critical coupling r\ ~ 4 for the RBEC state is consis- 
tent with our analytical result 



Vc 



2(kj_ 
7r V m 



In ( A/m + \J (A/to) 2 + 1 



(68) 



derived in |l9|. The difference between NBEC and RBEC 
states is that the chemical potential fi is of the order of 
to in the NBEC region but approaches zero in the RBEC 
region. 

The critical temperature, plotted as the solid line in 
Figj^k-, shows significant change from the weak to strong 
coupling. To compare it with the standard critical tem- 
perature for the idea boson gas, we solve the eauation[22j 



d 3 q 

(2tt) 3 



(eq - Mb) 



-b\ 



Mb)] 



n B 



with 



(69) 

= ^/q 2 + , boson mass tob = 2/i and bo- 
son number 7Jb = n/2, and show the obtained critical 
temperature as dashed line in Figj2ji. In the weak cou- 
pling region T c is very small and agrees with the BCS 
theory. In the NBEC region T c changes smoothly and 
there is no remarkable difference between the solid and 
dashed lines. Around the coupling rj c — 4, T c increases 
rapidly and then varies smoothly again. In the RBEC 
region, the critical temperature from our calculation de- 
viates significantly from the standard critical tempera- 
ture for ideal boson gas. Note that, T c is of the order 
of the Fermi kinetic energy e f ~ kj/(2m) in the NBEC 
region but becomes as large as the Fermi energy Ef in 
the RBEC region. The pseudogap A pg at T — T c , shown 
in Figl^t, behaves similarly as the critical temperature. 
To see clearly the pseudogap region, we present in Fig(2k 
the limit temperature T* as a dotted line. The pseudo- 
gap exists between the solid and dotted lines and begin 
to vanish on the dotted line. 

To explain why the critical temperature in the RBEC 
region deviates remarkably from the standard one for 
ideal boson gas, we calculate the boson number frac- 
tion ?'b = n&/(n/2) and the fermion number fraction 
rp = 1 — tb and show them as functions of the coupling rj 
in Fig|3] While there are only bosons at T c in the NBEC 




FIG. 2: The critical temperature T c (a), chemical potential 
/x(T c ) (b) and pseudogap A pg (r c ) (c) as functions of coupling 
r\ at A/m = 10 and kf/m = 0.5. T c , ^ and A pg are all scaled 
by the Fermi energy Ef. The dashed line is the standard 
critical temperature for the ideal boson gas in (a) and stands 
for the position y, = m in (b), and the dotted line in (a) 
is the limit temperature T* where the pseudogap starts to 
disappear. 



region, t-q is obviously less than 1 in the RBEC region. 
This conclusion is consistent with the results from the 
NSR theorypj, 03 ■ In the NBEC region, the binding 
energy of the molecules is Ef, ~ 1/too 2 = 2rj 2 £f, which 
is much larger than the critical temperature T c ~ 0.2e/, 
and the molecules can be safely regarded as point bosons 
at temperature near T c . However, the critical tempera- 
ture in the RBEC region is as large as the Fermi energy 
Ef, which is of the order of the molecule binding energy 
Eb — 2m. Due to the competition between the condensa- 
tion and dissociation of composite bosons in hot medium, 
the molecules can not be regarded as point bosons and 
the critical temperature should deviates from the result 
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FIG. 3: The boson number fraction re and the fermion num- 
ber fraction tf at the critical temperature T c as functions of 
the coupling 77 at A/m = 10 and kf/m = 0.5. 



for ideal boson gas. This may be a general characteris- 
tic of a composite boson system, when the condensation 
temperature T c is of the order of the molecule binding 
energy. The phenomenon can be explained by the com- 
petition between free energy and entropy jl7|: In terms 
of entropy a two-fermion state is more favorable than a 
one-boson state, but in terms of free energy it is less fa- 
vorable. Since the condensation temperature T c in the 
RBEC region is of the order of (n^/m^) 1 / 2 ~ (n/fi) 1 / 2 , 
we conclude that only for a system with sufficiently small 
value of kf/m, the standard RBEC critical temperature 
can be reached and is much smaller than 2m. 



V. APPLICATION TO MASSLESS FERMIONS: 
COLOR SUPERCONDUCTIVITY 

As a natural application of the relativistic GqG for- 
malism, we calculate in this section the transition tem- 
perature and pseudogap in two flavor color superconduc- 
tivity at moderate baryon density. The two flavor color 
superconducting quark matter corresponds to the ultra- 
relativistic case with n 3> m 3 , where m is the current 
quark mass. At moderate baryon density, the quark en- 
ergy gap due to color superconductivity is of the order of 
100 MeV, which is not located in the weak coupling re- 
gion. As a result, the pseudogap effect is expected to be 
significantly important near the critical temperature. To 
apply the present theory directly, we employ the gener- 
alized Nambu-Jona-Lasinio(NJL) model with scalar di- 
quark channel, which has been widely used to study color 
superconductivity at moderate baryon density. The La- 
grangian density is defined as 

C = $(j-fd M -m)il> + G B [(^) 2 + (^ 7 5^) 2 ](70) 
+G d (^7 5 T 2 A a C^ T ) (V T « 7 5 r 2 A a ^) , 

a=2,5,7 

where ip an d V' denote the quark fields with two flavors 
(Nf — 2) and three colors (N c = 3), Ti(i = 1, 2, 3) are the 



Pauli matrices in flavor space and A a (a = 1,2, ...,8) are 
the Gell-Mann matrices in color space, and G s and Gd 
are coupling constants in meson and diquark channels. 

At moderate baryon density, the chiral symmetry has 
already been restored and we need not consider the chiral 
condensate (ipip)- Since the current quark mass m is 
about 5 MeV, the quarks are nearly massless. The order 
parameter field for color superconductivity is defined as 



$ Q = -2G d ^ T C7z 7 5 r 2 A a V. 



(71) 



To simplify the calculation, one usually considers a spon- 
taneous color breaking from the SU(3) symmetry to a 
SU(2) subgroup. Due to the residual color SU(2) sym- 
metry, the effective potential in mean field approximation 
should depend only on the combination A 2 + A5 + A 2 
with A a = (3>a), and we can choose a specific gauge 
A sc = A 2 7^ 0, A5 = A7 = without loss of generality. 
In this gauge, the red and green quarks participate in the 
condensation, but the blue one does not. 

The detailed formalism of the GoG theory in the NJL 
model is similar to what we shown in sections [TT] and 
HIT] but becomes somewhat complicated due to the pres- 
ence of color and flavor degrees of freedom. Comparing 
the quark propagator in the NJL model with the one 
shown in above sections, the relativistic GoG scheme can 
be directly applied to the study of color superconduc- 
tivity, provided that we consider carefully the difference 
between the pairing including a blue quark and the pair- 
ing with only red and green quarks. The dispersion for 
red and green quarks is identical with the one obtained in 
the toy model, their excitation gap is A = (A 2 c +A 2 g y/ 2 , 
and the pair susceptibility x(l) should be multiplied by a 
factor Nf(N c — 1) where Nf and N c are flavor and color 
numbers of quarks. The new feature is that a gapless 
blue quark in the naive BCS mean field theory obtains a 
gap A pg in the GoG scheme. This can be understood by 
the fact that, the color symmetry is controlled only by 
the order parameters themselves, and fluctuations of any 
order parameter field $ a do not change it. At and above 
the critical temperature, A sc = 0, the color symmetry 
is restored, all colors become degenerate, and their gaps 
are just the pseudogap. 

The two flavor quark matter may exist in the region 
of fi = 350 — 500 MeV, where the strange quarks are not 
yet excited. Unlike the study in above sections in the 
canonical ensemble with fixed fermion number, people 
usually investigate color superconductivity in the grand 
canonical ensemble with fixed quark chemical potential. 
In this case, the quark number is not directly coupled to 
the calculation of the order parameter A sc and pseudogap 
Apg. For numerical calculations, we take the current 
quark mass m — 5 MeV, the often used quark momentum 
cutoff A = 650 MeV, and a fixed quark chemical potential 
fi = 400 MeV. We have checked that a reasonable change 
in the value of fj, does not bring qualitative difference. As 
is conventionally considered in the literatures, we use the 
pairing gap Ao at zero temperature to reflect the strength 
of the diquark coupling constant Gd- 
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In Fig[4]we show the critical temperature T c as a func- 
tion of A in the GqG theory and in the BCS mean field 
theory. While the critical temperature is not strongly 
modified by the diquark fluctuations in a wide range of 
Ao, the difference between the two is up to 20% in the 
strong coupling region with Ao — 200 MeV. In FigEl we 
show the pseudogap A pg at the critical temperature T c . 
In a wide range of the coupling, the pseudogap is of the 
order of 100 MeV, which is as large as the diquark con- 
densate A sc at zero temperature. Such a behavior means 
that the two flavor color superconductivity at moderate 
density is in the BCS-BEC crossover region and quite like 
the high temperature superconductivity in cuprates[l,0]- 
Since A sc vanishes at T = T c , the large pseudogap will 
bring significant effect at and above T c , such as the non- 
Fermi liquid behavior. In FigjSJ we show the temperature 
dependence of the diquark condensate A sc and pesudo- 
gap Apg at two values of Ao. With increasing tempera- 
ture, while the diquark condensate decreases, the pseudo- 
gap increases from zero. At low temperature, especially 
at zero temperature, we can safely neglect the pseudogap. 
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FIG. 4: The phase transition temperature T c for two flavor 
color superconductivity as a function of the diquark conden- 
sate Ao at zero temperature in the BCS mean field theory 
(dashed line) and in the GqG theory with diquark fluctua- 
tions (solid line). 



FIG. 6: The diquark condensate A sc (dashed lines) and pseu- 
dogap Apg (solid lines) in two flavor color superconductivity 
as functions of temperature scaled by T c for Ao = 100 MeV 
(upper panel) and 200 MeV (lower panel). 
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FIG. 7: The temperature dependence of the pseudogap A pg 
scaled by Aj /4 T 3/4 for A = 100 MeV and 200 MeV. 
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FIG. 5: The pseudogap A pg in two flavor color superconduc- 
tivity at the critical temperature T c as a function of Ao. 



While the pseudogap is small at low temperature and 
dominates the system only near and above T c , the di- 
quark fluctuations bring significant contribution to ther- 
modynamics at any temperature. In the low tempera- 
ture region, the temperature behavior of the pseudogap 
is significantly important, since it can tell us whether the 
coefficient Z\ or Zi dominates the pair fluctuations. In 
FigJT] we show the pseudogap at low temperature. In 
the region of T/T c < 0.1, it obeys a perfect power law 
A pg oc T 3 / 4 , which means that Z\ is the dominant one 
for the pair susceptibility. 

Considering the uncondensed diquarks, the total ther- 



11 



modynamic potential f2 can be expressed as 

^ ^cond ~l~ ^quark ~t~ ^diquark? (72) 

where the condensate and quark contributions take the 
same form as in the BCS theory, and the diquark contri- 
bution can be written as 

^diquark = M 1 - 4G dX(<z)]- (73) 

Since the coefficient Z\ controls the pair fluctuations at 
low temperature, the specific heat satisfies the power 
law C cx T 3 / 2 . As we mentioned above, the diquark 
contribution can be neglected only at sufficiently weak 
coupling. While the color superconductor at moderate 
baryon density may not reach the BEC condition, the ef- 
fect of diquark fluctuations on the thermodynamics may 
be remarkable, and it may bring significant astrophysi- 
cal consequences, such as the cooling process in compact 
stars. 



VI. SUMMARY 

We have generalized the non-relativistic GqG for- 
malism of BCS-BEC crossover to rclativistic fermion 
systems. The theory can describe the superfluid- 
ity/superconductivity with arbitrary strength of attrac- 
tive interaction, both in the symmetric phase and sym- 
metry breaking phase. The beyond-BCS effect at strong 
coupling brings in thermally excited bosons and con- 
tributes a pseudogap to fermion excitations. In such a 



formalism, we confirmed that there exists a BCS-NBEC- 
RBEC crossover in relativistic fermion systems. 

For color superconductivity at moderate baryon den- 
sity, while the BEC state can not be reached, the effect of 
diquark fluctuations is still remarkable and the naive BCS 
mean field theory breaks down when the temperature is 
close to the critical value. We investigate the two fla- 
vor color superconductivity at quark chemical potential 
/i = 350 — 500 MeV where the gap at zero temperature is 
of the order of 100 MeV. We found that the beyond-BCS 
effect strongly suppresses the transition temperature, and 
the pseudogap is very large near the critical temperature. 
This may strongly modify the thermodynamics of quark 
matter and bring significant astrophysical consequences 
in the study of compact stars. 

Such a theory can be applied to not only diquark con- 
densate ((qq)) at finite baryon density but also chiral con- 
densate ( (qq) ) at finite temperature and pion superfluid- 
ity at finite isospin density. The observation of oo bound 
states in strongly coupled quark-gluon plasma[26[ and a 
large thermal quark mass above the chiral phase tran- 
sition temperature in lattice QCD[27| indicate strongly 
the significance of the qq Bose-Einstein condensation and 
the quark pseudogap effe ct Il2l . [l3| . The study in this di- 
rection is under progress [28J. 

Acknowledgement: We thank Dr.Mei Huang for 
drawing our attention to [6j and [7J] and Gaofeng Sun 
for a careful check on the calculations in the Appendix. 
This work is supported by the grants No.NSFC10428510, 
No. 10575058 and the major state basic research develop- 
ing program 2007CB8 15000. 



APPENDIX A: PAIR SUSCEPTIBILITY AND ITS EXPANSION COEFFICIENTS 

In this appendix, we evaluate the explicit expression of the pair susceptibility and its momentum expansion. Com- 
pleting the trace in Dirac space and the Matsubara summation over the fermion frequencies, we obtain from the 
equations ([30]), (16]) and dTTJ) 



x(q) = 



d 3 k 

(2tt) 3 

/ 1 



' 1 - f(E k ) /Kq-k) £ k + £ k ) - /(£q-k) £ k "d/l k • (q - k) 



/(£k)-/(C- k )g k 

± t± 



k ' - /(Cq-k) #k + £k 



k ^k ~ S q -k + * 



2EZ 



2ekC q -k 



k-(q-k) 



2e k e 



q-k 



(Al) 



Taking its first and second order derivatives with respect to qo, we have 



Zi = 

z 2 = 



l-/(£ k )-/(£k) , f(E k )-.f(^) 



EZ + £ 



d 3 k 1 
(2?r) 3 2E- 
d 3 k 1 



f(EZ)-f(^) 



(E u 



Ck) 2 



(e£,£^ez,$), 



(A2) 
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Using the relation (E^ ) 2 — (Ck) 2 = ^ 2 > the coefficients can be rewritten as 



Zi 
Z 2 



1 



2A 2 
1 



n-2 
d 3 k 



d 3 k 



2A 4 J (2tt 



(/(£)-/(£)) 

TO 2 + «r v> 



(2tt) 3 



+ ( 



(A3) 



Taking the second order derivative of the susceptibility x with respect to q, we obtain another coefficient 



d 3 k 

(2tt) 3 



2E: 



i n-/(^k)-/(C) /(£ k )-/(£ k ) 



(^ 



C) 2 



C) 2 



2r(e k : 

A 2 
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l-2/(££)\ e 



37 +« 



2^r 



37 -G 



2Sr 



2K 



2e 4 

Zfc k 



(A4) 



with a; = cos# and f'(x) being the first order derivative of the Fermi-Dirac distribution. 

In the non-relativistic limit with kf <C to, \fj, — m\, A <C m, all the terms including anti-fermion dispersions can be 
safely neglected, and the relativistic dispersions are reduced to 6c = k 2 / (2m) — (p — m) and E k = + A 2 . Taking 
into account |q| <C to, we have 



Xnr(<?) 



d 3 k 

(X) 3 



1 - /(gk) - /gq-k) 3c + 6c /(3c) ~ /(gq-Q 3 k ~ gk 

3 k + Cq-k - qa 2E k E k - 6i-k + go 2^ 



(A5) 



which is just the same as the one given in [1,0], and the expansion coefficients Zi, Zi and £ 2 are reduced to@, @] 



Z 2 = 
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In the RBEC limit with \i — ► 0, we can expand Zi in powers of /i, Zi ~ R[i + 0(/i 3 ), with the expansion coefficient 
R given by 



R = 



rf 3 k 
(2T) 3 
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where Ek = \/e£ + A 2 is the dispersion at [i = 0, and Zi and £ 2 can be simplified 



as 



Z, = 



d 3 k 



A 4 7 (2tt) 3 
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